Mathematical Biosciences and Engineering

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has continuously evolved since its emergence in the human population in 2019. As of 1st August 2025, more than 1,700 Omicron subvariants have been designated by the Pango nomenclature system. The Pango nomenclature system designates a new lineage based on genetic and epidemiological information of SARS-CoV-2 strains. However, there is a possibility that strains that have similar genetic backgrounds and the same phenotype are given different Pango lineage names. In this paper, we propose a new algorithm, called FindPart-w, which can identify groups of viral lineages that share the same relative effective reproduction numbers. We introduced a new lineage replacement model, called the constrained RelRe model, which constrains groups of lineages to have the same relative effective reproduction numbers. The FindPart-w algorithm searches the equality constraints that minimise the Akaike Information Criterion of constrained RelRe models. Using hypothetical observation count data created by simulation, we found that the FindPart-w algorithm can identify groups of lineages having the same relative effective reproduction number in a practical computational time. Applying FindPart-w to actual real-world data of time-stamped lineage counts from the United States, we found that the Pango lineage nomenclature system may have given different lineage names to SARS-CoV-2 strains even if they have the same relative effective reproduction number and similar genetic backgrounds. In conclusion, this study showed that viruses that had the same relative effective reproduction number were identifiable from temporal count data of viral sequences. These findings will contribute to the future development of lineage designation systems that consider both genetic backgrounds and transmissibilities of lineages.

Human immunodeficiency virus (HIV) persistence remains a major barrier to cure due to the existence of long-lived latent reservoirs that evade immune clearance and persist despite combination antiretroviral therapy (ART). Although ART effectively suppresses viral replication, treatment interruption often leads to rapid viral rebound originating from these latent reservoirs. In this study, we develop a deterministic mathematical model describing the in vivo dynamics of HIV infection incorporating uninfected CD4+ T cells, infected cells, latent reservoirs, deep latent reservoirs, and infectious and non-infectious virions, while explicitly accounting for the therapeutic effects of reverse transcriptase inhibitors (RTIs), protease inhibitors (PIs), and Tat transcription inhibitors. Analytical results establish positivity and boundedness of solutions and derive the effective reproduction number Re using the next-generation matrix approach. Stability analysis shows that the virus-free equilibrium is locally asymptotically stable when Re < 1, while viral persistence occurs when Re > 1. Numerical simulations were performed to investigate therapy interactions, viral rebound following treatment interruption, and the impact of drug efficacy on viral set-points and latent reservoir dynamics. To further explore therapy interactions, three-dimensional viral set-point surfaces and heat maps were generated to examine how combinations of infection inhibition, viral production inhibition, and transcriptional inhibition influence viral dynamics. The simulations reveal that Tat inhibition suppresses viral transcription, thereby reducing the transition of infected cells into productive infection and limiting viral propagation when combined with conventional ART mechanisms. The therapy parameter planes further demonstrate that strong transcriptional inhibition promotes the transition of infected cells into deep latency, supporting the emerging block-and-lock strategy for functional HIV cure. In addition, a three-dimensional eradication boundary surface and therapy cube were constructed to identify regions of parameter space where Re < 1, corresponding to successful viral control. These visualizations show that viral eradication is unlikely when therapies act independently but becomes achievable when multiple therapeutic mechanisms act simultaneously. Overall, the results highlight the critical role of transcriptional inhibition through Tat-targeting therapies in complementing existing ART regimens. By simultaneously suppressing viral replication and promoting deep latency, Tat-based combination strategies may significantly reduce viral rebound and contribute to long-term functional control of HIV infection.

Against the backdrop of global climate change and accelerating population mobility in 2025, chikungunya fever (CHIKF) exhibited a trend of worldwide spread, significantly increasing the difficulty of controlling tropical mosquito-borne diseases. To enhance the precision of intervention strategies, this study developed an age- and sex-structured human-mosquito interaction dynamic model based on data from the largest CHIKF outbreak ever recorded in China, and conducted a targeted analysis of prevention and control strategies. By decomposing the basic reproduction number and examining population heterogeneity, asymptomatic males aged 15-59 years were identified as the core transmission group. Optimal control analysis revealed that the synergistic implementation of three measures-- reducing the effective human-to-mosquito transmission rate, reducing the effective mosquito-to-human transmission rate, and suppressing mosquito population density--could reduce the overall infection rate by 95.7586%. Among these, mosquito population suppression should be prioritized as a universal core strategy; however, its protective effect on females aged 60 years and above was relatively weak, warranting particular attention. The study further demonstrated that asymmetric intensity combinations targeting these three intervention pathways--such as intensity profiles of "10%, 90%, 90%" or "60%, 80%, 90%"--could achieve effective outbreak control. This research elucidates population-specific transmission patterns and key pathways for intervention intensity, providing a theoretical and strategic foundation for the precise control of mosquito-borne diseases. It also provides actionable operational insights to support rapid response and strategy optimization for future emerging outbreaks. Author summaryCHIKF is a mosquito-borne viral disease that is gradually spreading from tropical regions to other areas. To achieve more precise control of this disease, we developed an age- and sex-structured analytical model based on the largest CHIKF outbreak in China, aiming to provide a scientific basis for responding to potential future outbreaks with inherent uncertainties. The study found that asymptomatic males aged 15-59 years were the primary drivers of transmission and should be prioritized as a key population for reducing viral spread in prevention efforts. When evaluating the effectiveness of different intervention strategies, females aged 60 years and above were the least affected by the implemented measures, indicating that this group should strengthen personal protection to lower their infection risk. Among all control measures, mosquito suppression was the most effective, suggesting that vector control strategies should be prioritized in future outbreak responses.

Within many bacterial colonies, persister cells exist as a subpopulation that is tolerant to antibiotics and other stressors, yet not genetically distinct from the rest of the colony. A recent study has proposed epigenetic inheritance as a mechanism that leads to the presence of persister cells. We analyze a nonlocal PDE--ODE model introduced in that study to describe the epigenetic inheritance process and establish its mathematical well-posedness, including existence, uniqueness, and nonnegativity of solutions. We identify a sharp parameter threshold delineating extinction from persistence of the colony: below this threshold the washout equilibrium is globally asymptotically stable, while above it a unique positive equilibrium exists and the population is weakly persistent. Notably, this threshold is independent of the internal community structure.

Variability in gene expression among single cells and growing cell populations can arise from the stochastic nature of protein synthesis, which often occurs in random bursts. This study investigates the variability in the expression of a growth-sustaining protein, whose concentration is regulated by a negative feedback loop due to cell growth-induced dilution. We model the distribution of protein concentration using a Chapman-Kolmogorov equation for single cells and a population balance equation for growing cell populations. For single cells, we derive an explicit solution for the protein concentration distribution in state space and represent it as a Bessel function in Laplace space. For growing populations, we find that the distribution satisfies a Heun differential equation with singular boundary conditions. By addressing the central connection problem for the Heun equation, we quantify the population-level protein distribution and determine the Mathusian parameter, which characterizes population growth. This work provides a comprehensive analytical framework for understanding how stochastic protein synthesis impacts gene expression variability and population dynamics.

Mass-conserved reaction-diffusion systems are used as mathematical models for various phenomena such as cell polarity. Numerical simulations of this system present transient dynamics in which multiple stripe patterns converge to spatially monotonic patterns. Previous studies indicated that the transient dynamics are driven by a mass conservation law and by variations in the amount of substance contained in each pattern, which we refer to as "pattern flux". However, it is challenging to mathematically investigate these pattern dynamics. In this study, we introduce a reaction-diffusion compartment model to investigate the pattern dynamics in view of the conservation law and the pattern flux. This model is defined on multiple intervals (compartments), and diffusive couplings are imposed on each boundary of the compartments. Corresponding to the transient dynamics in the original system, we consider the dynamics around stripe patterns in the compartment model. We derive ordinary differential equations describing the pattern dynamics of the compartment model and analyze the existence and stability of equilibria for the reduced ODE with respect to the boundary parameters. For a specific parameter setting, we obtained results consistent with previous studies. Moreover, we present that the stripe patterns in the compartment model are potentially stabilized by changing the parameter, which is not observed in the original system. We expect that the methodology developed in this paper is extendable to various directions, such as membrane-induced pattern control.

Infectious disease surveillance systems in tropical countries show that respiratory disease incidence generally manifests as year-round activity with weak fluctuations and irregular seasonality. Previously, using a ten-year time series of influenza-like illness (ILI) collected from outpatient clinics in Ho Chi Minh City (HCMC), Vietnam, we found a combination of nonannual and annual signals driving these dynamics, but with unknown mechanisms. In this study, we use seven stochastic dynamical models incorporating humidity, temperature, and school term to investigate plausible mechanisms behind these annual and nonannual incidence trends. We use iterated filtering to fit the models and evaluate the models by comparing how well they replicate the combination of annual and nonannual signals. We find that a model including specific humidity, temperature, and school term best fits our observed data from HCMC and partially reproduces the irregular seasonality. The estimated effects from specific humidity and temperature on transmission are nonlinearly negative but weak. School dismissal is associated with decreased transmission, but also with low magnitude. Under these weak external drivers, we hypothesize that stochasticity makes a strong sub-annual cycle more likely to be observed in ILI disease dynamics. Our study shows a possible mechanism for respiratory disease dynamics in the tropics. When the external drivers are weak, the seasonality of respiratory disease dynamics is prone to the influence of stochasticity.

Mainland China experienced multiple waves of COVID19 pandemic during 2020 2022, driven by emerging variants and changes in public health and social measures (PHSMs). We developed a hypergraph-based Susceptible Vaccinated Exposed Infectious Recovered Susceptible (SVEIRS) model to reconstruct epidemic dynamics across 31 provinces, capturing transmission heterogeneity associated with clustered contacts. We assessed key characteristics of transmission at national and provincial levels during four outbreak periods: initial, localized predelta, Delta, and widespread Omicron, which accounted for 96.7% of all infections. We found significant diversity in transmission contributions across cluster sizes, with a small fraction of larger clusters responsible for a disproportionate share of infections. Counterfactual analyses showed that reducing clustersize heterogeneity, while holding overall exposure constant, could have lowered national infections by 11.70 to 30.79%, with the largest effects during Omicron period. Ascertainment rates increased over time but remained spatially heterogeneous with a range: (14.40, 71.93)%. Population susceptibility declined following mass vaccination (to 42.49% in Aug 2021, nationally) and rebounded (to 89.89% in Nov 2022) due to waning immunity with variations across the provinces. Effective reproduction numbers displayed marked temporal and spatial variability, with higher estimates during Omicron. Overall, these results highlight critical role of group contact heterogeneity in shaping epidemic dynamics.

Background Six years after its emergence, SARS-CoV-2 continues to have a substantial burden. The impact of vaccination and the optimal timing of its rollout remain uncertain given existing population immunity and variability in outbreak timing between summer and winter. Methods The US Scenario Modeling Hub convened its 19th round of ensemble projections for COVID-19 hospitalizations and deaths in the United States, where eight teams projected trajectories in each US state and nationally from April 2025 to April 2026 under five scenarios regarding vaccine recommendations and timing. Recommendations had two eligibility scenarios (high-risk individuals only and all-eligible) and two timing scenarios (classic start: mid-August, earlier start: late June). These were crossed to create four scenarios and were compared against a counterfactual scenario with no vaccination. Findings Compared to no vaccination, our ensemble projections estimated 90,000 (95% PI 53,000-126,000) hospitalizations averted in the high-risk and classic timing scenario across the US. Expanding to all-eligible age-groups averted an additional 26,000 (95% PI 14,000-39,000) hospitalizations, which when coupled with the early vaccination timing, was projected to further reduce national hospitalizations by 15,000 (95% PI -3,000-33,000). The majority of teams projected both summer and winter waves. Implications We project COVID-19 will cause significant hospitalizations and deaths in the US in the 2025-26 season and estimate significant benefits from a broad all-eligible vaccination recommendation. The results also suggest an additional benefit is likely to be gained from an earlier vaccination campaign. Funding Centers for Disease Control and Prevention; National Institute of Health (US), National Science Foundation (US)

Understanding the dynamics of HIV epidemics is important to control them effectively. Classical methods that mainly rely on occurrence data are limited by the fact that an unknown part of the epidemic eludes sampling. Since the early 2000s, phylodynamic methods have enabled the estimation of key epidemiological parameters from virus genetic sequence data. These methods have the advantage of being less sensitive to partial sampling and to provide insights about epidemic history that even predates the first samples. In this study, we analysed 2,205 HIV sequences from the French ANRS PRIMO C06 cohort. We identified and were able to reconstruct the temporal dynamics of two large clades that represent the HIV-1 epidemics in the country. Using Bayesian phylodynamic inference models, we found that the first clade, from subtype B, originated in the end of 1970s, grew rapidly during the 80s before decreasing from 2000 to 2015 and stagnating since then. The second clade, from circulating recombinant form CRF02_AG, emerged and spread in the 80s, grew again in the early 2000s, before declining slightly. We also estimated key epidemiological parameters associated with each clade. Finally, using numerical simulations, we investigated prospective scenarios and assessed the possibility to meet the 2030 UNAIDS targets. This is one of the rare studies to analyse the HIV epidemic in France using molecular epidemiology methods. It highlights the value of routine HIV sequence data for studying past epidemic trends or designing public health policies.

Plasmodium vivax transmission from humans to mosquitoes depends on the density of gametocytes that in turn depends on asexual parasite replication and gametocyte commitment. These processes are often analyzed separately, despite being biologically linked and measured with substantial uncertainty. We used a joint Bayesian latent-variable model to simultaneously analyze parasite density, gametocyte density, and mosquito infectivity while accounting for measurement error and propagating uncertainty across linked processes. The model was applied to individual-level data from three P. vivax transmission studies conducted in Ethiopia (n = 455). A tenfold increase in gametocyte density was associated with more than a twofold increase in the odds of mosquito infection (odds ratio [OR] = 2.32, 95% credible interval [CrI]: 2.12-2.54). Asexual parasite density was also positively associated with infectivity after accounting for gametocyte density (OR = 1.74, 95% CrI: 1.60-1.90), and inclusion of parasite density improved predictive performance. Pathway decomposition within the joint model indicated that approximately 41% of the parasite-infectivity association operated through gametocyte density. Increasing age was associated with lower asexual parasite density but higher gametocyte density, resulting in minimal overall association with infectivity. Predicted infection probability increased sigmoidally with gametocyte density, remaining low at lower densities before increasing sharply and approaching a plateau at higher densities. Gametocyte density produced the largest predicted changes in the proportion of infected mosquitoes, while asexual parasite density added predictive information not fully captured by measured gametocyte density alone. This approach links molecular parasite measurements with mosquito infection risk while accounting for measurement uncertainty and provides an interpretable framework for studying the P. vivax infectious reservoir. Author SummaryMalaria transmission occurs when mosquitoes ingest sexual-stage parasites, called gametocytes, during a blood meal. In Plasmodium vivax infections, human-to-mosquito transmission depends on linked biological stages, including asexual parasite replication, gametocyte production, and mosquito infection. These processes are closely connected and often measured with uncertainty, making them difficult to study using standard approaches that analyze them separately. In this study, we applied a joint Bayesian model that analyzes parasite density, gametocyte density, and mosquito infectivity together while accounting for uncertainty in laboratory measurements. Using data from three studies in Ethiopia, we quantified how parasite density, gametocyte density, and host characteristics relate to mosquito infection. The analysis showed that measured gametocyte density alone did not fully explain variation in infectivity, and that asexual parasite density provided additional predictive information. We also found that age was associated differently with asexual parasite and gametocyte densities, resulting in little overall association with infectivity. This approach helps link molecular parasite measurements with transmission outcomes and improves understanding of the P. vivax infectious reservoir in endemic settings.

Oscillatory behaviour is important in multiple biological contexts. However, the inherent nonlinearity and high dimensionality of mathematical models in biology makes proving the existence and the localization of limit cycle oscillations challenging. Here, we provided an elementary proof for the existence and a method for rigorously localizing the oscillatory solutions in a class of benchmark biomolecular oscillators. To construct the proof, we used a geometric approach based on Brouwers Fixed Point theorem. We constructed a toroidal-like manifold within a positively invariant set by removing the hypervolume containing the fixed point and the trajectories converging to it along its stable manifold. We showed that the vector field describing the system dynamics maps a cross section of the toroidal-like manifold onto itself. The existence of a limit cycle solution in this manifold was guaranteed by Brouwers Fixed Point theorem. For different sets of initial conditions in these cross-sections, we used an interval-based Reachability Analysis to localize the oscillatory behaviour that complements the Brouwers Fixed Point theorem approach. These results add a simple and elegant approach to demonstrating the existence of limit cycles in biomolecular systems as well as a method for rigorous localization of the region of existence.

In "Experiment-free exoskeleton assistance via learning in simulation", Luo et al. [1] present an ambitious framework for developing exoskeleton controllers through reinforcement learning exclusively in computer simulation. The authors report that a control policy trained on a small dataset from one subject was directly transferred to physical hardware, reducing human metabolic cost during walking, running, and stair climbing by more than any prior device. If confirmed, this would represent a major breakthrough for the field of wearable robotics and their clinical applications. However, a close examination of the published materials casts doubt on these claims. The reported experimental results violate physiological limits on the relationship between mechanical power and muscle energy use during gait2,3,4. The algorithmic claims are surprising and cannot be verified; in contrast with established replicability standards in machine learning5,6, executable code has not been made available. We conclude that the goals of this study have not yet been verifiably achieved and make recommendations for avoiding publication errors of this type in the future.

Cerebral blood flow (CBF) control is essential for normal brain function and is disrupted in pathological conditions. Arterial diameters are tightly regulated to provide on demand increases in blood flow in regions of neuronal activity. Pericytes (PCs) exhibit robust myogenic tone and may also respond to neuronal activity to fine-tune local resistance and blood flow. Thus, mural control of microcirculatory resistance may extend beyond arteries and arterioles. Yet, PCs electrophysiology and contractility have not been thoroughly characterized, and this prohibits an integrated view of brain blood flow control. In this study, we develop a detailed mathematical model of mural cell electrophysiology, Ca2+ dynamics and biomechanics. The model is informed by electrophysiological data in smooth muscle cells (SMCs) or PCs and predictions are compared against pressure-induced responses in isolated arterioles and capillaries, respectively. Simulations recapitulate myogenic constrictions and examine differences in contractile dynamics as we move from arterioles to proximal and distal capillaries. In arteriole-to-capillary transitional (ACT) zone PCs, increased mechanosensitivity, more Ca2+ influx through non-selective cation (NSC) channels and/or a higher sensitivity of the contractile apparatus to Ca2+ can compensate for reduced L-type voltage-operated (VOCC) Ca2+ influx and allow for robust constrictions at the lower operating pressures of capillaries relative to the arterioles. A significant Ca2+ influx through NSC relative to VOCC, however, can decouple the PCs contractile apparatus from electrical signaling. Vasoactivity to chemomechanical stimuli along the arteriole to capillary axis is progressively driven by VOCC-independent Ca2+ influx and Ca2+ sensitization with slow kinetics. The proposed cell model can form the basis for detailed multiscale and multicellular models that will examine physiological function at a single vessel or vascular network levels and investigate CBF control in health and in disease. Key pointsO_LIA mural cell model of electrophysiology, calcium (Ca2+) dynamics and biomechanics is informed by data and adapted for modeling cerebral arteriole smooth muscle cells and capillary pericytes. C_LIO_LIIon channel activities are characterized by patch-clamp electrophysiology in isolated cerebral smooth muscle cell and pericytes, and capillary and arteriole electromechanical responses to transmural pressure changes are assessed using novel ex vivo preparations. C_LIO_LIMyogenic constrictions in arterioles can be reproduced by pressure-induced non-selective cation channel (NSC) activation that depolarizes the cell, opens L-type Ca2+ channels (VOCCs) and increases Ca2+ influx. C_LIO_LIRobust myogenic constrictions in arteriole-to-capillary transition (ACT) zone pericytes may reflect significant Ca2+ influx through NSC, increased mechanosensitivity, or higher sensitivity of the contractile apparatus to Ca2+, potentially compensating for reduced VOCC density relative to arteriolar smooth muscle. C_LIO_LIA significant contribution of NSC relative to VOCC in Ca2+ influx, can decouple the contractile apparatus from electrical signaling. C_LIO_LIThe model shows how gradients in ionic activities, mechanosensitivity and/or Ca2+ sensitivity can alter contractile phenotype and electromechanical coupling along the arteriole to capillary continuum. C_LIO_LIThe proposed model can form the basis for detailed multiscale and multicellular models that will investigate cerebral blood flow control in health and in disease. C_LI

Pulmonary embolism (PE) is a sudden blockage of lung arteries, usually caused by a blood clot that travels from the deep veins of the legs. As the world becomes more sedentary and lifestyle diseases emerge, deaths from PE are expected to rise in the next 20 years. For instance, the United States records annual deaths of 60 per 100,000 people. The degree to which these deaths are affected by demographic, socioeconomic and environmental predisposing factors as well as how they vary across time and space remains an open science question. In this paper, we conduct a detailed statistical and spatial-temporal study PE mortality counts across US counties from 2005 to 2022. Our study shows that study shows that PE mortality is not randomly distributed in space and time but concentrated in most counties in Arkansas, Mississippi, Kansas, Missouri, Oklahoma, Louisiana, Nebraska, Tennessee, and Texas. We also established that age is a statistically significant predictor (mean coefficient of 0.52) of PE mortality especially in counties of Mississippi, Kansas, Missouri, Tennessee, Illinois, Kentucky, Texas and Virginia. Our results thus provide empirical support for prioritizing regionally targeted PE prevention policies. Furthermore, the adopted county-level analysis uncovered granular geographic patterns that are usually obscured in state or national level analysis. Our study thus provides actionable evidence to support geographically tailored strategies aimed at reducing mortality by pinpointing counties with consistently elevated PE mortality risk at different timescales.

The role of passive muscle mechanical properties in explosive performance is important to understand to maximize training and performance in strength sports. The objective of this study was to create and test a Structural Equation Model (SEM) to analyze the interrelations between muscle mechanical properties and neuromuscular performance in competitive weightlifters. The MyotonPRO was used to measure muscle stiffness, muscle tone, muscle elasticity, relaxation time and creep of four major muscles: Quadriceps Femoris, Hamstrings, Trapezius and Biceps Brachii, on thirty elite male weightlifters. This involved performance metrics of the rate of force development (RFD), countermovement jump (CMJ) and time to contraction threshold (TCT). AMOS was used to analyze the direct and indirect relationships between variables through SEM analysis. These findings indicated that muscle stiffness and relaxation time were significant predictors of explosive performance measures (p < 0.05) but there were weak or no relationships between tone, elasticity and creep. The model proposed had good fit indices (CFI = 0.97, RMSEA = 0.049) indicating its structural soundness. These results present the significance of muscle stiffness and relaxation time as important predictors of neuromuscular performance. The proposed model suggests an effective structure of monitoring athletes, their performance optimization, and individual training design in strength-based sports.

Entomological surveillance plays a crucial role in areas where malaria remains endemic, yet gathering data on mosquito populations is often expensive and complicated, particularly in remote locations with challenging logistics and inconsistent sampling schedules. Access to extensive time series data on mosquito species at specific sites would greatly enhance insights into seasonal trends and the biting habits of vectors of malaria parasites. Gaps in mosquito count records pose a significant challenge for researchers and public health officials seeking to establish early warning systems and effective vector control programs. In this study, we apply quantitative machine learning techniques to address missing data in estimates of mosquito abundance collected from 2009 to 2016 in Bolivar State, Venezuela. We evaluated Linear Regression, Stochastic Linear Regression, K Nearest-Neighbor, and Gradient Boosting methods for imputing missing counts of Anopheles mosquitoes, employing a leave-one-out cross-validation strategy. Additionally, we developed a predictive malaria transmission model incorporating mosquito abundance and climate variables (El Nino 3.4 Index, rainfall, and mean air temperature) as covariates. Our generalized time series model forecasts malaria incidence of Plasmodium vivax and Plasmodium falciparum based on climate dynamics and imputed mosquito data. Model performance was assessed using root mean square error, mean absolute error, and mean absolute percentage error. The final results demonstrated that machine learning imputation significantly improved the accuracy and reliability of P. vivax malaria incidence predictions but failed to predict P. falciparum incidence. The study demonstrates that method choice significantly influences the reconstruction of seasonal abundance patterns and the performance of malaria incidence models. Nevertheless, the proposed models strengthen the foundation for targeted interventions and surveillance in endemic regions. Despite limitations in data continuity and coverage, the findings highlight the value of combining multiyear entomological data sets with robust imputation and sensitivity analyses to improve predictive modeling in resource-constrained, malaria-endemic settings.

In October 2025, mpox virus clade I infections have been detected among men who have sex with men (MSM) in the EU/EEA, suggesting local transmission in MSM sexual networks. Given the large outbreak of mpox among MSM in 2022 and the uncertain transmission parameters of clade I in the European context, clade I poses a public health concern to the EU/EEA. This work assesses the potential effect of increasing the mpox vaccine uptake among MSM via two contributions. First, building on the European MSM and Trans Persons Internet Survey 2024, we estimate the mpox vaccine uptake among MSM as well as the proportion who are unvaccinated but willing to get vaccinated for 28 countries in the EU/EEA. Specifically, we fit Bayesian mixed-effects models for the vaccine and recovery status of an individual depending on their number of sexual partners and country. Second, we develop a susceptible-infectious-recovered model on a sexual contact network to estimate the reduction of the reproduction number if vaccines are provided to MSM who are willing to get vaccinated. Our results suggest a substantial willingness for mpox vaccination among MSM if mpox cases increase and a large reduction of the effective reproduction number if this willingness is met. These findings highlight a large potential of increasing mpox vaccine uptake among MSM and preventing future mpox outbreaks in the EU/EEA.

The basic reproduction number of an infectious disease is known to depend on the structure of contacts between individuals in a population. This relationship has been explored mathematically through two well-known models: one which depends on a matrix of contact rates between different demographic groups, and another which depends on the variability of contact rates over the population. Here we introduce a model that combines and generalises these two approaches. We derive a formula for the basic reproduction number and validate it through comparisons to simulated outbreaks. Applying this method to contact survey data collected in Belgium between 2020 and 2022, we find that our model produces higher estimates of the basic reproduction number and larger relative changes over periods when social contact behaviour was changing during the COVID-19 pandemic. Our analysis suggests some practical considerations when using contact data in models of infectious disease transmission.

Cancer-associated fibroblasts (CAFs) are a central component of the tumor microenvironment that facilitate a supportive niche for cancer progression and metastasis. Experimental evidence suggests that CAFs can facilitate estrogen-independent tumor growth, thereby reducing the efficacy of anti-hormonal therapies. Understanding and quantifying the complex interactions between tumor cells, hormonal signalling, and the microenvironment are crucial for designing more effective and individualized treatment strategies. We propose a mathematical framework to explore the influence of CAFs on ER+ breast cancer progression and to evaluate strategies to mitigate their impact. We develop a system of nonlinear ordinary differential equations that substantiates the experimental observations by providing a mechanistic basis for the role of CAFs in regulating estrogen-independent growth dynamics. We then employ optimal control theory to evaluate distinct therapeutic approaches involving monotherapy or combinations of: (i) inhibition of tumor-to-CAF signaling, (ii) inhibition of CAF-to-tumor proliferative signaling, and (iii) endocrine therapy. Taken together, our results demonstrate that CAF-targeted strategies can enhance treatment efficacy across various estrogen dosing regimens. Our study provides new insights into the potential of CAF as a therapeutic target that could help to improve existing approaches for endocrine therapies.